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We study two-particle long-range rapidity correlations arising in the early stages of heavy ion 
collisions in the saturation/Color Glass Condensate framework, assuming for simplicity that one 
colliding nucleus is much larger than the other. We calculate the two-gluon production cross section 
while including all-order saturation effects in the heavy nucleus with the lowest-order rescattering 
in the lighter nucleus. We find four types of correlations in the two-giuon production cross section: 
(i) geometric correlations, (ii) HBT correlations accompanied by a back-to-back maximum, (iii) 
away-side correlations, and (iv) near-side azimuthal correlations which are long-range in rapidity. 
The geometric correlations (i) are due to the fact that nucleons are correlated by simply being 
confined within the same nucleus and may lead to long-range rapidity correlations for the produced 
particles without strong azimuthal angle dependence. Somewhat surprisingly, long-range rapidity 
correlations (iii) and (iv) have exactly the same amplitudes along with azimuthal and rapidity shapes: 
one centered around A<f> = n with the other one centered around A<j> = (here A<j> is the azimuthal 
, angle between the two produced gluons). We thus observe that the early-time CGC dynamics in 

D ■ nucleus-nucleus collisions generates azimuthal non-flow correlations which are qualitatively different 

Ph ' from jet correlations by being long-range in rapidity. If strong enough, they have the potential of 

mimicking the elliptic (and higher-order even-harmonic) flow in the di-hadron correlators: one may 
need to take them into account in the experimental determination of the flow observables. 
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I. INTRODUCTION 



^ PACS numbers: 25.75.-q, 25.75.Gz, 12.38.Bx, 12.38.Cy 

43 . 

Long-range rapidity correlations between pairs of hadrons produced at small azimuthal angles with respect to each 
' other were discovered recently in heavy ion (AA) [H-Hl], proton-proton (pp) [f|, and proton-nucleus (pA) collisions 
, Q . Due to the particular shape of the corresponding correlation function, with a narrow correlation in the azimuthal 
■ angle A<f> and a wide correlation in pseudo-rapidity separation At/, these correlations are often referred to as the 
0\ ; "ridge". 

There appears to be a consensus in the community that the origin of these long-range rapidity correlations is in 
the very early-time dynamics immediately following the collision. A simple causality argument demonstrates that a 
correlation between two hadrons produced far apart in rapidity may arise only in their common causal past, that is, 
in the early stages of the collision [7, 8]. However, the detailed dynamical origin of these "ridge" correlations is not 
, completely clear. 

It has been proposed in the literature [7H"l5| that the "ridge" correlations may arise in the classical gluon field 
dynamics of the parton saturation physics/Color Glass Condensate (CGC). (For reviews of saturation/CGC physics 
see [1614201 ].) Indeed classical gluon fields, which in the McLerran-Venugopalan (MV) model [2l| - [23j domina te g luon 
production in heavy ion collisions, do lead to a rapidity-independent distribution of the produced gluons [24 - |27l | over 
C3 • rapidity intervals of up to Ay < l/a s , which is the upper limit of their validity (with a s the strong coupling constant). 
Correlations between such classical fields, introduced in the process of averaging over their color sources, do have a 
long range in rapidity 0, [H, [Hj]. Moreover, it was observed in 0, [TlHIH that the diagrams giving rise to such rapidity 
correlations also lead to a narrow correlation in the azimuthal direction, in qualitative agreement with the shape of 
the "ridge" correlation. 

One has to keep in mind that in heavy ion collisions the early-stage azimuthal correlation may be washed out by 
the final state interactions leading to thermalization of the produced medium and its hydrodynamic evolution. (Note 
though, that the rapidity correlation is not likely to be strongly affected by such late-time dynamics.) It was argued, 
however, that the effect of the radial flow in the hydrodynamic evolution of the quark-gluon plasma (QGP) would be 
to (re-)introduce the azimuthal correlations @. 

Another potential complication with the CGC explanation of the "ridge" is the fact that rapidity-dependent correc- 
tions to classical gluon fields do become important at rapidity of the order of Ay ~ l/a s . These corrections come in 
through the non-linear Balitsky-Kovchegov (BK) |29l - l32| and Jalilian-Marian-Iancu-McLerran-Weigert-Leonidov- 
Kovner (JIMWLK) (33l-[36| evolution equations. Rapidity-dependent corrections are very important for describing the 
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hadron multiplicity distribution in rapidity, dN/dy, in the CGC framework [371 ] . As dN/dy docs depend on rapidity 
rather strongly in RHIC heavy ion data, and also strongly (albeit less so) in the LHC data, it is very hard to describe 
without the rapidity-dependent nonlinear evolution. It is possible that similar rapidity-dependent corrections may 
significantly affect and potentially destroy the long-range structure of the rapidity correlations due to classical gluon 
fields of the MV model. Note that progress on this issue has been made in [1JJ, indicating that inclusion of small- a; 
evolution still leaves the "ridge" reasonably flat in rapidity until Ay ~ l/a a when the correlation disappears. 



To elucidate the above questions and concerns, and to improve the precision of the CGC predictions for the "ridge" - 
like correlations, it is important to be able to calculate the two-particle correlation in the CGC framework beyond the 
lowest order. While some works do consider the role of small- a; evolution and multiple rescatterings in the correlation 
function [111, [l3|, [38l - l40l | , most of the phenomenological approaches [TJ, [HI, l4fl - l43T ] simply include the saturation effects 
into the lowest-order calculation by evolving the unintcgrated gluon distributions with the running-coupling BK 
(rcBK) nonlinear evolution |44l - l46j |. 



The aim of this paper is to begin to analytically include saturation effects into the two-gluon correlation function 
in nucleus-nucleus collisions. Indeed full inclusion analytic of saturation effects originating in both nuclei would be a 
very hard problem: even the sin gle gluo n production in the quasi-classical MV limit of AA collisions can be dealt with 
only numerically at present [26|, [23, E3, |48[ . To make the problem more tractable we assume that one of the colliding 
nuclei is much larger than the other o ne, such that saturation effects are important only in interactions with the larger 
nucleus. In a more formal language (49|, for the quasi-classical MV model, we assume that 1 <C A\ -C A2, where 
A\ and A2 are the atomic numbers of the smaller and larger nuclei, while o? s A 1 ^ 3 ~ 1 and a 2 s A\' 3 <C 1. Equivalently 
one can say that we are interested in production of gluons with transverse momentum kx > Q s i, where Q s \ is the 
saturation scale of the smaller nucleus. We do not impose any constraints on Ht compared to the saturation scale Q S 2 
of the larger nucleus (Q S 2 ^ Qsi ^ Aqcd)- We stress here that this setup is not what is usually referred to as the 
pA collision in the saturation/CGC terminology: since A\ 3> 1, the two gluon production cross section is dominated 
by the gluons produced in interactions of different nucleons in the smaller nucleus, which may also be viewed as a 
lowest-order saturation correction. 



Below we calculate the two-gluon production cross section and the corresponding correlation function including the 
all-order saturation effects in the larger nucleus, while keeping them to the lowest non-trivial order in the smaller 
nucleus. The paper is structured as follows. In Sec.|IT]we set up the problem. In the process of properly defining the 
correlation function we unexpectedly find a correlation originating in a simple fact that the nucleons in both nuclei are 
confined to within the nuclear radii. The effect is quite generic, and is not specific to the saturation/CGC approach 
used in the calculations in much of this paper: this correlation should be present in any model of AA collisions which 
properly takes into account the geometry of the collision. We refer to these correlations as the "geometric correlations" . 
We argue that this effect gives a non-trivial contribution to the correlator even if the two-gluon production cross section 
is given by disconnected Feynman diagrams, i.e., the correlation is not generated dynamically in the interactions. This 
correlation may be long-range in rapidity, though it does not have any non-trivial azimuthal structure. 



The two-gluon production cross section in a heavy-light ion collision is calculated in Sec. Mil in terms of the Wilson 
line correlators. We show that the two-gluon production cross section is related to various correlators of four adjoint 
Wilson lines: the adjoint dipole, quadrupole, and a double-trace operator. These correlators contain both the multiple 
rescatterings of the quasi-classical MV approximation, along with the small- x BK/JIMWLK evolution. 



We study the long-range rapidity correlations in Sec. I IV! We evaluate the Wilson-line correlators describing the 
interaction with the target nucleus using the quasi-classical MV approximation and the large- N c limit. By analyzing 
the two-gluon correlation function at the lowest non-trivial order, and in agreement with the earlier calculations in 
the literature [TJ EH HH we nn d both the near-side (A</> = 0) and away-side (A</> = it) long-range rapidity 
correlations. (Here A<p is the azimuthal angle between the momenta of the produced gluons.) The two correlations 
are identical in their azimuthal shapes (as functions of Acf>), such that the correlator can be expanded into a Fourier 
series in terms of only the even harmonics cos 2 n Acf>. We discuss the possibility that such non-flow correlation may 
complicate experimental extraction of the contribution of the true QGP flow to the flow observables V2n- In addition to 
the geometric and long-range rapidity correlations, the obtained two-gluon production cross section contains Hanbury- 
Brown-Twiss (HBT) correlations [51| ■ A somewhat peculiar feature of this HBT correlation is that it is accompanied 
by an identical back-to-back peak as well. We conclude in Sec. [V] by stressing the difference between near-side 
and away-side correlations calculated here and the mini-jet correlations: the near-side correlations in this work are 
long-range in rapidity, while the jet near-side correlations are local in rapidity. 
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II. CORRELATION FUNCTION AND "GEOMETRIC" CORRELATIONS 



Definition of the correlator 



Following a standard approach used in experimental analyses of particle correlations [H, 52] the correlation function 
can be defined as 



dN 12 



c(k uyi ,k 2 ,y 2 )=Ar d J N dyid2k2 d T 1 (1) 



d' 2 k±dyi d 2 k-2dy2 



where 

dN 1 da 



d 2 kidyi a ine i d 2 k 1 dyi 

and 

dN 12 1 da 



d 2 kidyid 2 k 2 dy 2 a lnel d 2 k 1 dyi d 2 k 2 dy 2 



(2) 



(3) 



are the single- and double-particle multiplicity distributions with ai ne i the net inelastic nucleus-nucleus scattering 
cross section. Two-dimensional transverse vectors are denoted by v = (v x ,v y ) with their length vt = \v\. The 
normalization factor Af in Eq. ([1]) is fixed by requiring that the number of particle pairs measured in the same 
("real") event N 12 is equal to the number of ("mixed") pairs with particles coming from different events (TV) 2 . Here 
we are interested in Arj-A(j) correlations with A</> = 4>\ — <j) 2 the difference between the azimuthal angles <f>i, <p 2 of the 
two momenta of the produced particles, and 6rj = r\\ — r/ 2 ~ y\ — y% the difference between the pseudo-rapidities 771, r\ 2 
of the two particles. For such a correlation, with the magnitudes of the transverse momenta k\T and k 2 T constrained 
to some chosen data bins, the normalization factor is fixed by 

N [ dfa dyi # 2 dy 2 -5— i§? A = f Widyi „, , [ #2 dy 2 , 2 f , ■ (4) 
J d 2 k 1 dyid 2 k 2 dy 2 J d*k x dyi J d z k 2 dy 2 

(We assume for simplicity that the number of produced particles is very large, N S> 1, such that N — 1 « TV and a ine [ 
is the same in both Eqs. @ and © since the production cross section of producing exactly one particle is negligible.) 
Combining Eqs. ([T]), ©, ([3]), and ((4]) we rewrite the correlation function in terms of cross sections as 



I^xdyx^-Jd^dy,^ 



da 

d' 2 k\dyi d 2 k2<!>j2 



c(k uyi ,k 2 ,y 2 ) = ^ r - — r - 1 - ( 5 ) 



/ d(f>x d Vl dcj) 2 dy 2 dikldyi % k2dy2 



da da 
d' 2 k±dyi d 2 k-2dy2 



As outlined in the Introduction, we take one of the nuclei (the projectile) to be much smaller than the other, 
Ai <C A 2 , such that saturation effects in it are minimal and nucleons in this ( "first" ) nucleus interact with the other 
("second") nucleus independent of each other. At the same time, the projectile nucleus is still large enough, A\ S> 1, 
such that the two-gluon production is dominated by the gluons generated by the interactions of two different nucleons 
from the first (projectile) nucleus with the second nucleus. (This interaction with the two nucleons in the first nucleus 
can be considered a lowest-order saturation effect.) Denote the 1- and 2-nuclcon wave functions of the projectile 
nucleus A\ by 'fj(b) and ^ n(b\,b 2 ): they are normalized such that 

J d 2 b\^f I (b)\ 2 = A 1 , J d 2 b 1 d 2 b 2 \y II (b u b 2 )\ 2 = A 1 (A 1 -l)^A 2 . (6) 

With the help of these wave functions the single- and double-gluon production cross sections can be written as 

r d 2 Bd%\* I (B-b)\ 2 (J^\ (7a) 



d 2 kdy J \d 2 kdyd 2 b 

= I d 2 Bd 2 bl d% 2 |*„(B — &!,£?— b 2 )\ 2 I J U T^ -J^^ 



d 2 k\dy\d 2 k 2 dy 2 J \d 2 kidy\d 2 bi d 2 k 2 dy 2 d 2 b 2 



where B is the impact parameter between the two nuclei and 6, 61 , b 2 the transverse positions of the nucleons in the 
projectile nucleus, all measured with respect to the center of the second (target) nucleus, as shown in Fig. Q] for the 
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FIG. 1. Transverse plane geometry of the two-particle production in the collision of a smaller projectile nucleus (^4i) with a 
larger target nucleus (A2). The two smaller circles represent two nucleons in the nucleus A\ (see text for details). 



two-gluon production process. (Transverse vector b labels the position of the incoming nucleon in the single gluon 
production case, while b\ and hi label positions of projectile nucleons for two-gluon production.) Here 

da pA2 



dPkdydPb 



is the cross section for the gluon production (with fixed transverse momentum fc, rapidity y, and transverse position 
b) in the collision of a nucleon (p) with the target nucleus. The angle brackets (. . .) in (JT]) denote averaging in the 
target nucleus wave function along with summation over all the nucleons in the target nucleus [2l| - [23l . I33I436I . |49| . 
Eqs. ([7]) can be used in Eq. ([5]) to give us the two-gluon correlation function in the heavy-light ion collision 



C(k 1 ,yx,k 2 ,y 2 ) = 



J d 2 Bd 2 b x dfadyi \Vi(B -bx)\' 



vrf 2 fei dyi d 2 b x 



-) fd 2 Bd 2 b 2 d<p 2 dy 2 \y I (B-b 2 )f 



1 da pA 2 

\ d 2 k2 dy2 d 2 b2 ' 



/ d 2 B d 2 h d 2 b 2 d(j>i dyx d<t> 2 dy 2 \^fu(B — bi,B — b 2 )\ 2 
Jd 2 Bd 2 b 1 d 2 b 2 \V n (B - 61,-B- b 2 )\ 2 



dijP A 2 d gf A2 
d 2 k\ dyi d' 2 bi d 2 k2 



da pA2 dcrPA2 . > 
d 2 k\dy\d' 2 b\ d 2 k2dy2d 2 b2 



J d 2 Bd 2 b 1 \Vi(B - 61) 



d<yP A 2 
d 2 k\ dyi d 2 bi 



Jd 2 Bd 2 b 2 \^i(B-b 2 )\ 2 



d<yP A2 
d 2 k2 dy2 d 2 b'2 



(8) 



Here we model the nucleus as a ba g of independent nucleons, which is a correct description at the leading order in 
the atomic number A [21M23I , l49l IHol. |53| . In such case the single- nucleon light-cone wave function squared is simply 
equal to the nuclear profile function in the projectile nucleus []J 

\9 I (b)\ 2 =T 1 (b). (9) 

(The nuclear profile function for a nucleus with density p(b, z) is defined by the integral over the longitudinal coordinate 



OO 

T(b)= J dzp(b,z), 



(10) 



such that for a spherical nucleus of radius R and constant density p it is T(b) = 2 p v ' R 2 — b 2 . 
Without any loss of generality we can perform the integral over B in Eq. (|7a[) obtaining 



da 



d?k dy 



= A x / d 2 b 



daP A2 



d 2 kdyd 2 b 



(11) 



1 We do not show the spin and isospin indices explicitly in the wave functions: in our notation the wave function squared is implicitly 
averaged over all nucleon polarizations, since both colliding nuclei are unpolarized. 
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For a sufficiently large projectile nucleus, A\ ^> 1, one can assume that the two-nucleon wave function can be 
factorized, 



®ii(bi,b 2 ) = $r(6i)*/(6 2 ), 
such that, with the help of Eq. (|9]) we can write 

|*//(bi,b 2 )| 2 =Ti(&i)Ti(&a). 

Using this in Eq. (ffb| we get 

f d 2 Bd 2 h d 2 b 2 T 1 (B - b 1 )T x {B - b 2 ) 



d 2 k\dy\ d 2 k 2 dy 2 
Substituting Eqs. dTTJ) and ([H]) into Eq. © yields 



C{k u yi 1 k 2l y 2 ) 



d 2 k\dy\d 2 b\ d 2 k 2 dy 2 d 2 b 2 



4 2 d 2/2 ( d i k % y2 2 d 'i b2 ) 



Jd 2 B d 2 b Y d 2 6 2 d<fc dyi # 2 dy 2 Ti(B - b^T^B - b 2) v d . fc2 dy2 d . &2 



da pA 2 



da pA 2 



Jd 2 B d 2 h d 2 b 2 Ti(B -b 1 )T 1 (B -b 2 ) 



d 2 kidy\d' 2 bi d? '■k^dy^d^bi 



(12) 

(13) 
(14) 



(15) 



To complete the calculation one needs the single- and double-gluon production cross sections which, when used in 
Eq. (fTB"j) , would give us the correlation function. Before we proceed to construct them, let us study a simple example 
elucidating the nature of one of the correlation types contained in correlator (fT5j) . 



B. Geometric correlations 



Let us consider the simplest possible example of particle (gluon) production mechanism where the interaction of 
the two nucleons in the first nucleus with the second nucleus in Eq. (|14[) factorizes, 



da pM 



da pA2 



d 2 k\dy\d 2 b\ d 2 k 2 dy 2 d 2 b-2 



da pA2 



d 2 k\dy\d 2 b\ I \d 2 k 2 dy 2 d 2 b 2 



da pA2 



(16) 



This contribution comes from the disconnected Feynman diagrams and is usually identified as the uncorrelated part 
of the two-gluon production cross section. However, it is clear that substituting Eq. (fT6"| into Eq. (|15[) does not reduce 
the correlation function to zero: instead one gets 



C(ki,yi, k 2 ,y 2 ) = 



fd 2 h 



dyi 



' da pA 2 \ 
^d' 2 ki dyi d' J bi I 



Jd 2 b 2 



4 2 d 2/2 ( d i k d 2 a d P y 2 2 d 2 b2 ) 



J d 2 B d 2 b x d 2 b 2 #i dy x dcj> 2 dy 2 T x {B - 6 X ) T x (B - b 2 
J d 2 B d 2 h d 2 b 2 Ti(B -b 1 )T 1 (B-b 2 ) 



da pA 2 
WkTd~yTcFb 1 

dcr pA2 
d 2 k\dy\d 2 b\ 



da pA 2 
d' 2 k 2 dy2d 2 b2 



da pA 2 
d 2 k-2dy 2 d 2 b2 



Jd 2 h 



dcrP A 2 
d 2 k\ dy\ d 2 b\ 



Jd 2 b 2 



da pA 2 

d 2 k 2 dy 2 d 2 b 2 



(17) 



which, in general, could be non-zero. 

Certainly if the 6-dependence factorizes from the rapidity and azimuthal dependence in the cross section 

da pA2 



d 2 kdy d 2 b 



(18) 



then the correlation function (JTTJ) is zero: however, such factorization is not always the case. For gluon production 
in the saturation framework, the cross section is a complicated function of kT/Q s (b,y), which means it is not in a 
factorized form and thus the correlator (|17[) is not zero. Note that in the MV model (which does not contain the small- 
x evolution), gluon production is rapidity- independent, and, if one neglects the dependence of the gluon production 
cross section on the angle between k and b, the correlator (|17p becomes zero. (Dependence of gluon production cross 
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section on the collision geometry in the MV approximation is not very strong, peaking at non-pcrturbatively low 
momenta [54|.) 

For the general case in Eq. (|17j) we observe a possible non-trivial correlation in the two-gluon production described 
by disconnected Feynman diagrams. If the gluon production cross section (|18[) is a slowly varying (but not constant) 
function of rapidity, as is the case in the saturation/CGC framework near mid-rapidity, this correlation would be long- 
range in rapidity. The origin of this correlation is somewhat peculiar: even though the two-nucleon wave function in 
Eq. (|13p is factorized and, hence, represents uncorrelated nucleons, these two nucleons are correlated by the simple 
fact of being parts of the same bound state, the projectile nucleus. In other words, the probability of finding two 
nucleons at the impact parameters b\ and b 2 is proportional to 

~ J d 2 BT 1 (B-b 1 )T 1 (B-b 2 ) (19) 

and is not a product of two independent probabilities after all impact parameters B of the incoming nucleus are 
integrated over: this is a correlation. Note also that the presence of real wave-function correlations, that is, non- 
factorizable corrections to the right-hand-side of Eq. (fill)) , would also lead to some nontrivial two-particle correlations 
in Eq. (JT7J). 

If we define the correlation function at the fixed nuclear impact parameter B by not integrating over B in Eqs. ([7]) 
and using the result in Eq. (|5|), we get 



C(k 1 ,y 1 ,k 2 ,y 2 ;B) = 



J d%, dft dyi \mb M 2 U?CW" 




J d 2 b 2 dfo dy 2 \MB - b 2 )| 2 U?Ofa>" 




/ <Ph d% 2 deb, d Vl dfa dy 2 \* n (B -b u B- b 2 )\* { d 4%% bi d ^^J_ 





J dH, d% 2 \9 n (B -b^B-b^ 

X / d^lvKfl-bOl 8 { d -4t\W ) Id 2 b 2 \^(B-b 2 )\^ ( a^pg ) 

(20) 

One can see that this fixed-impact parameter correlation function in Eq. ([20]) is zero, C(B) = 0, for the factorized 
wave function from Eq. (|13j) and for disconnected-diagram interactions from Eq. (|16|) . Thus di-gluon correlations 
due to our "geometric" correlation mechanism seem to also disappear when the impact parameter is fixed exactly. 
However, such precise determination of the impact parameter is impossible in an experimental analysis, where one is 
able to fix the collision centrality \B\ in a certain interval, but one can not fix the direction of B. The integration 
over any range of the impact parameter \B\ (or the integration over the angles of B keeping \B\ constant) is likely to 
introduce these geometric correlations, as follows from Eq. (|15l) . Note that the presence of non-trivial wave function 
correlations, i.e., correlations beyond the factorization approximation in Eq. (|13p . may also lead to correlations which 
may survive in Eq. (|20[) even for a fixed impact parameter B and uncorrelated interactions (|16j) . 

Despite its simplicity, the non-vanishing correlation in Eq. (|17p is one of the main results of this work. In the 
saturation/CGC framework it may lead to long-range rapidity correlations similar to the observed "ridge" correlation. 
Indeed azimuthal correlations are missing in Eq. (j 1 T[) : such correlations may be formed in heavy ion collisions due to 
radial flow, as was argued in Q. Indeed a lot more work is needed to compare this result with experiment. 



III. TWO-GLUON PRODUCTION AND CORRELATIONS 



In this Section we are going to calculate the two-gluon production cross section for the heavy-light ion collisions 
working in the saturation/CGC framework. As described above we assume that A 2 ^> A\ 3> 1, such that the 
saturation effects are resummcd to all orders only in the target nucleus with atomic number A 2 . While the saturation 
effects are not very important in the projectile nucleus with the atomic number A\ , the fact that A\ 3> 1 implies that 
the two gluons arc predominantly produced in collisions of different nucleons in the projectile nucleus with the target 
nucleus. 

According to the standard technique (HI l20j , the calculation will proceed using the light-cone perturbation theory 
(LCPT) [55l [56j to construct the gluon wave functions in the incoming projectile nucleus. The production cross 
section will be obtained by convoluting the square of this light-cone wave function with the interaction with the target 
nucleus, which will be described by Wilson lines along the light cone [Tt], [29J, which include the saturation effects in 
the form of both the Glauber-Mueller multiple rescatterings [53[ and the BK/JIMWLK evolution equations. 

The diagrams contributing to the scattering amplitude of the two-gluon production in the collision of two nucleons 
from the projectile nucleus with the target nucleus are shown in Fig. [5J There the two nucleons in the incoming nucleus 
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are denoted by horizontal solid lines representing two quarks in the wave functions of those nuclcons: throughout this 
work we will model the nucleons by a single valence quark each. (Generalization of our results to a more realistic 
nucleon wave function is straightforward.) Interaction with the target, which happens over a time-scale much shorter 
than the time scale of preparing gluons in the wave functions, is denoted by a vertical dashed line. 



bi 



b 2 




i 



D E F 

FIG. 2. Diagrams describing the two-gluon production in the heavy-light ion collision. Two horizontal solid lines denote valence 
quarks inside the two nucleons in the projectile nucleus. Vertical dashed line denotes the interaction with the target nucleus, 
while the vertical dotted lines denote intermediate states. 



A. Light cone wave functions of the gluons 



An essential part of the calculation is finding the light cone wave functions of the two interacting nuclcons generating 
gluons in the incoming nucleus. The diagrams in Fig. [2] include vertical dotted lines denoting the intermediate states 
contributing light-cone energy denominators to the wave functions. The energy denominators, and hence the wave 
functions of the two nucleons, factorize. While this is obvious in diagrams C and D in Fig. [2j we will explain this 
factorization for other diagrams in a little more detail below. 

We assume that the projectile nucleus is moving along the light cone "+" direction and the nucleons carry very 
large p + = (p° +p 3 )/v / 2 momenta. At the same time, the gluons carry light cone momenta kf and k£ such that 
p + S> k^ , fcj • The gluons will be produced with rapidities far away from the fragmentation regions of the colliding 
nuclei. The small light cone momenta of the gluons make them dominate the energy denominators, generating a much 



larger contribution than that of the (valence) quarks [57 1 



The only difference between the diagrams A and B in Fig. [2] is the ordering of the emitted gluo ns. This will 
only affect the energy denominators that result from the light-cone perturbation theory rules (55l l56l| . When these 
two diagrams are added together these energy denominators factorize. Denoting the light cone energies of the two 
gluons in diagrams A and B by E\ and E 2 such that E\ = k\ T j (2 k±) and E2 = k% T /(2 k£) we see that adding the 
contributions of the energy denominators in those graphs gives 

1 1 + 1 1 11 

Ei Ei + E2 E2 Ei + E2 Ei E2 

The result factorizes into the terms describing one or another gluon. Analysis of the diagrams E and F leads to the 
same factorization of energy denominators [58| . 

The two-gluon light cone wave function for diagrams A+B and E+F is 

V Qb (asi, b u A; x 2 , b 2 , A') = ^(xi, 61, A) i>\x 2 , b 2 , A') = - 9 - 2 (t a )i (t») 2 ^ = ^ \f x (a . 2 ~ ^ \f' (22) 

\xi-bi\ z \X2~b2Y 

where g is the QCD coupling constant and the polarization vector is e\ = — (l/\/2)(A, i) [55|,[56|. The factors of (t a )i 
are the SXJ(N C ) generators in the fundamental representation with N c the number of quark colors, while the subscript 
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1, 2 denotes the associated nucleons. The single nucleon's soft gluon wave function is [5j 

7r \x — or 



(23) 



The diagrams C+D have a wave function different from Eq. (|22[) by a minus sign. The minus sign difference will be 
absorbed in the interaction term. 



B. Single gluon production 



Our calculation for the two-gluon production cross section will be closely following that for the single gluon pro- 
duction in a collision of a single nucleon with a nucleus (a pA collision). For this reason, and also because we need 
the single inclusive gluon production cross section to construct the correlator (|15|) . let us briefly review the results. 

The singlc-gluon inclusive production cross section in pA collisions is illustrated in Fig. [3J which shows the square 
of the scattering amplitude. The cross denotes the measured gluon. The vertical solid straight line in Fig. [3J is the 
final-state cut, while the vertical dashed lines denote interactions with the target just like in Fig. [2] 







FIG. 3. Diagrams contributing to the square of the scattering amplitude for the single gluon production in pA collisions. The 
cross denotes the measured produced gluon. 



The gluon production cross section can be written as [59 62] (with the projectile proton represented by a single 
valence quark) 



(24) 



with the coordinate notation defined in Fig. [3j The wave function squared summed over polarizations and colors is 

^ ri , Mr . {yM ^^.^ 2 (25) 

with Cf = — l)/(2N c ) the fundamental-representation Casimir operator, while the interaction with the target is 



Inti{x,y,b) = 



Tr[U x Ul\ - Tr[U x ul] Tr[CW+] + 1^ 



AT? - 1 



(26) 



The latter object is written in terms of expectation values (in the target nucleus wave function) of the adjoint Wilson 
lines along the x + light cone: 



Pexp<z# J dx^ A~ {x + , x~ = 0, x) 



(27) 



where is the gluon field of the target in the adjoint representation. 
Defining the adjoint (gluon) dipole S'-matrix by 



Sg{xi,x 2 ,v) 



(28) 
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we rewrite Eq. ([24]) as 

(J^pi) = / M y e " <fe ' (x " a) ' if^ [Sg{x ^ v) - SG{x ' b ' y) ~ SG{b ^ v) + 1] • (29) 

Using Eq. (|2"9"|) in Eq. (fTTT) we obtain the gluon production cross section in the heavy-light ion collision. 

The dipole amplitude Sg in Eq. (|29|) contains the Glauber-Mueller multiple rescatterings [HH (recovered by putting 
y = in the arg ument of Sg), along with the energy dependence included though the BK/JIMWLK evolution 
equations [58l |62| : the cross section formula remains the same in both cases. The former case would correspond 
to consistently treating the problem in the MV model [59h61| . The latter case would include small- a; evolution 
corrections in the rapidity interval between the produced gluon and the target nucleus: one may worry that such 
inclusion would be asymmetric, since the evolution in the rapidity range between the projectile and the produced 
target is not included. A more symmetric treatment of the gluon production problem in pA collisions with evolution 
corrections included in all rapidity intervals is presented in (5o | (see also [l6l |20| for a more pedagogical presentation). 
Inclusion of evolution in the rapidity interval between the projectile and the produced gluon(s) is beyond the scope 
of the present work: it should be possible along the lines of [58| though. 



C. Two-gluon production with long-range rapidity correlations: "square" of the single gluon production 

The two-gluon production in heavy-light ion collisions is easily constructed by analogy to the single-gluon production 
calculation of Sec. IIII B[ The diagrams contributing to the square of the scattering amplitude for the double gluon 
production in heavy-light ion collisions arc shown in Fig. 01 written as a direct product of the gluon production 
processes in the interactions of each of the nucleons from the projectile nucleus with the target. Just as in Fig. [3] 
the vertical dashed lines in Fig. @] represent interactions with the target, while the vertical solid line denotes the final 
state cut. 



bi 




FIG. 4. Diagrams contributing to the two-gluon production cross section in the heavy-light ion collision. For clarity the 
diagrams are shown as a direct product of gluon production processes in collisions of the two interacting nucleons from the 
projectile nucleus with the target nucleus. 

The evaluation of the diagrams in Fig. 0] is straightforward. The two-gluon wave function squared is obtained by 
using Eq. (j25|) twice, which yields 

V' ab (^i,b 1 ,A; : E 2 ,6 2 ,A')V ab *(yi,b 1 ,A;y 2 ,b 2 ,A') 

A, A' ,a,6 

16 a 2 s C 2 F xi - bi yi - b x x 2 - b 2 y 2 - b 2 , . 

tt2 l^-bil 2 ' 12/! - foil 2 \x 2 -b 2 \ 2 ' |y 2 -b 2 | 2 [M) 
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with all the coordinate labels explained in Fig. |4j 

Interactions with the target are treated similarly. The traces of Wilson lines in the top and bottom parts of Fig. @] 
are exactly the same as in Eq. (|26| and appear to factorize, suggesting absence of dynamically generated correlations. 
This is not so. Note, that the averaging over the target is applied to both parts of Fig. 2] simultaneously: we thus get 
an averaged product of traces of Wilson lines, 

Int 2 (x 1 ,y 1 ,b 1] x 2l y 2l b 2 ) = ^j^-jTrp^U^] - j^-j Tr^^J - Tr^E/tj + A 

j^j Tr[U X2 Ul) - Tr[U X2 Ul) - j^-j Tr[U b2 Ul 2 ] + l)^ , (31) 

which does not, in general, factorize into the product of two separately target- averaged interactions from Eq. (l26l) . 
Using Eq. (|14p we write 

= r^m I d2B d% d 2 b 2 T i(B - b i) T dB - b 2 )d 2 x, d 2 Vl d 2 x 2 d 2 y 2 e"* "H'i-vi)-* fc-(- 3 -«i) 



da s 

d 2 k 1 dy 1 d 2 k 2 dy 2 ~ [2(27 

x ^ i' ab (x 1 ,b 1 ,X;x 2 ,b 2l X')^ ab *(y ll b 1 ,X;y 2 ,b 2 ,X')Int 2 (x 1 ,y ll b 1 ;x 2 ,y 2 ,b 2 ), (32) 



A, A' ,a,b 

which, with the help of Eqs. ([30]) and fl3U becomes (cf. jlol ]) 
da 



/ qu Z , = TTT I d 2 B d 2 h d 2 b 2 T^B b 1 ) T,{B b 2 ) d 2 x x d 2 Vl d 2 x 2 d 2 y 2 e ~< ^(—mM M» a -») 
a z k\dy\d* k 2 dy 2 lb7r s J 

xi - b\ yi - bi x 2 b 2 y 2 - b 2 



xi-bi\ 2 I2/1 — 61 1 2 \x 2 -b 2 \ 2 |i/2 
1 

i\L 2 - 



^-rr^C/tj _ _J_ r.^t/tj _ _L_ Tr[c/biC/ tj + ^ 

2^-7 ^[C/^C/tj - ^JL_ Tr[U X2 Ul 2 ] - j^-j Tr[U b2 Ul 2 ] + l\ V (33) 



This is the expression for the two-gluon production cross section coming from the diagrams in Fig. 2) The interaction 
with the target can be evaluated in the quasi-classical multiple rescattering approximation, as we will show later. The 
rapidity evolution can be included using the JIMWLK equation. Note, however, that when the rapidity difference 
between the two gluons is sufficiently large, \yi — y 2 \ > l/a s , one has to include the evolution corrections in the 
rapidity interval between the produced gluons, such that the Wilson lines in the two parenthesis in Eq. (|33[) should be 
taken at different rapidities. A similar effect had to be included in the two-gluon production cross section in DIS in 
[Hj]. In such regime one would also need to include evolution corrections in the rapidity window between the projectile 
and (at least one of) the produced gluons. Inclusion of the evolution corrections in terms of the weight functional W 
of the JIMWLK evolution equation into the two-gluon production cross section in nucleus-nucleus collisions was done 
in fill ]. In this work we will limit ourselves to the quasi-classical regime where no evolution corrections are required 
and the cross sections are rapidity-independent. 

Note that for scattering on a large nuclear target, if one resums powers of a 2 A 1 / 3 , as is the case in the Glauber- 
Mueller (GM) rescattcrings (and in the BK/ JIMWLK evolution for which the GM rescatterings serve as the initial 
condition), and if one takes the large- A c limit, the expectation values of the traces in Eq. (|33p factorize, such that 



(Trp^UlJ Tr[U X2 Ul 2 }) 



« (Tr[U Xl Ul}) {Tr[U X2 Ul 2 }) . (34) 

large- n c , largc-A 2 



This leads to factorization of Eq. (Tj"6")) being valid in the large-iV c and large-target-nucleus limit. As shown above, 
even in this factorized regime one may obtain a non-trivial correlation function due to the geometric correlations. 



D. Two-gluon production with long-range rapidity correlations: "crossed" diagrams 

Before we proceed to evaluating the correlations contained in the cross section (|33[) . let us point out another 
contribution to the two-gluon production cross section arising from squaring the sum of the diagrams in Fig. [2] When 
squaring the diagrams in Fig. [5] it is possible that the gluon emitted by one nucleon in the amplitude will be absorbed 
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by another nucleoli in the complex conjugate amplitude. The corresponding contributions to the cross section are 
shown in Figs. [5J and [5j where crossing gluon lines do not form a vertex. We will refer to these diagrams as the 
"crossed" graphs. The diagrams obtained from those in Figs. [5] and [6j by a mirror reflection with respect to the cut 
correspond to the 6i <-> b 2 interchange, and will be automatically included in the cross section to be calculated below 
since it will contain integrals over all b\ and b 2 . The sum over different orderings of gluon emissions in the (complex 
conjugate) amplitude of Fig. [6] is implied, but is not shown explicitly. 



bi 




FIG. 5. Diagrams contributing to the two-gluon production cross section, with the gluon emitted by each nucleon in the 
amplitude absorbed by another nucleon in the complex conjugate amplitude. The top cross denotes the gluon with momentum 
fci , while the bottom one denotes the gluon with momentum &2 . 



In evaluating the graphs in Figs. [5J and [5] we note that there are non-trivial color factors due to the interaction 
terms, which we absorb into the wave function. The expression for the wave function squared is 

16 a 2 s C F xt - bi y 2 - b 2 x 2 - b 2 yi - bi , . 

^ 2N C l^i — foil 2 ' |y 2 -b 2 | 2 \x 2 -b 2 \ 2 ' lyi-bil 2 ' 1 j 

Defining the color-quadrupole operator (63T - [67| 

Q( XuX2 , Xs , X4 ) = J J— (Tr[U Xl Ul 2 U x MS) (36) 

after some algebra the interactions with the target in both Fig. [5J and Fig. [6J can be shown to be equal to (with the 
terms ordered in the same way as the diagrams in Figs. [5J and [5] and the transverse coordinates defined in the upper 
left diagram of each of these two figures) 

Int crossed (x 1 ,y 1 ,b 1 ,x 2 ,y 2 , b 2 ) = Q(x 1 ,y 1 ,x 2 ,y 2 ) - Q(x 1 ,y 1 ,x 2 , b 2 ) - Q(x 1 ,y 1 ,b 2 ,y 2 ) + 5" G (a;i,yi) 

- Q(x 1: bi,x 2l y 2 ) + Q(x 1 ,b 1 ,x 2 ,b 2 ) + Q(x 1: b l7 b 2 ,y 2 ) - S G (xi, 61) 

- Q{t>i,yi,x 2 ,y 2 ) + Q{b 1 ,y 1 ,x 2 ,b 2 ) + Q(bi,yi,b 2 ,y 2 ) - SG(bi,yi) 

+ S G (x 2 ,y 2 ) - S G (x 2 ,b 2 ) - S G (b 2 ,y 2 ) + 1. (37) 

The only difference between the contributions of the diagrams in Figs. [5] and |6j is in the exponential factors for the 
Fourier transform into the transverse momentum space. 

The longitudinal momentum flow patterns in Figs. [5] and [6J are different from that in Fig. [4] This is illustrated in 
Fig. [71 which shows the flow of the "plus" momentum component through the first diagram in Fig. [51 Note that the 
change in the "plus" momentum component is negligible in the eikonal interactions with the target considered here. 
Requiring that the incoming quark lines carry the same "plus" momentum both in the amplitude and in the complex 
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FIG. 6. Another set of diagrams contributing to the two-ghion production cross section, with the gluon emitted by each nucleon 
in the amplitude absorbed by another nucleon in the complex conjugate amplitude. Again the top cross denotes the gluon with 
momentum fei, while the bottom one denotes the gluon with momentum k-2- Summation over the different orderings of gluon 
emissions (e.g. which gluon is emitted first or second) is implied but is not shown explicitly. 



conjugate amplitude, one would obtain = k% : however such requirement is not correct. The actual scattering 
happens between two nuclei, and it is the momenta of the whole incoming nuclei which have to be equal both in the 
amplitude and in the complex conjugate amplitude. Hence for ^ k^ the diagrams like that in Fig. [7] would only 
correspond to different redistributions of the projectile nucleus momentum between the nucleons in it in the amplitude 
and in the complex conjugate amplitude without changing the same "plus" momentum of the whole nucleus on both 
sides of the cut. Hence the k^ = k£ condition is not necessary for the "crossed" diagrams. 




Combining the factors in Eqs. (|35[) and (|37[) . and inserting Fourier transform exponentials, we arrive at the expression 
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d 2 k 1 dy 1 d 2 k 2 dy 2 [2(2? 



for the two-gluon production cross section contribution resulting from the "crossed" diagrams from Figs. [5] and [6] 

J d 2 B d 2 bi d 2 b 2 T\{B — bi) Ti(B — b 2 ) d 2 xi d 2 yi d 2 x 2 d 2 y 2 

16 a 2 C F x± - bi y 2 - b 2 x 2 - b 2 y 1 - b x 
, tt 2 2N C l^-bil 2 ' \y 2 -b 2 \ 2 \x 2 -b 2 \ 2 ' li/i-bil 2 

Q(x 1 ,y i ,x 2 ,y 2 ) - Q(x 1 ,y l7 x 2} b 2 ) - Q(x 1 ,y 1 ,b 2 ,y 2 ) + S G (x 1 ,y 1 ) - Q(x 1 ,b 1 ,x 2 ,y 2 ) + Q(x l7 b 1 ,x 2 ,b 2 ) 
Q{x 1 ,b 1 ,b 2 ,y 2 ) ~ S G {x 1 ,b 1 ) - Q(b 1 ,y 1 ,x 2 ,y 2 ) + Q(b 1 ,y 1 ,x 2 ,b 2 ) + Q(bi,y 1 ,b 2 ,y 2 ) - S G (bi,y{) + S G (x 2 ,y 2 ) 



e — i fei-(a)i— 1/2)— i fe2-(sc2— 3/1) _|_ g— i fei-(a!i— 2/2)+* fe2-(as2— yi) 



5 , G(^2,b2)-5 , G (b 2 ,y 2 ) + l 



(38) 



Just like in Eq. (j33j), the dipole and quadrupole scattering amplitudes in Eq. (|38j) can be evaluated either in the 
MV model or by using BK and JIMWLK evolution equations. The quadrupole amplitude evolution equation was 
derived in the large- N c limit in [63| . and beyond the large- iV c limit in j65|. Again, Eq. is valid only as long as 
the rapidities y\ and y 2 of the two produced gluons are close to each other, \y 2 — yi\ < l/a s , such that no small-x 
evolution corrections need to be included in the [yi , y 2 ] rapidity interval. 

E. Two-gluon production with long-range rapidity correlations: the net result 

Equations (|33p and (|38[) , when combined, give us the two-gluon production cross section in the heavy- light ion 
collisions: 

do d(j S q Uare d(T crosse( i (39) 



d 2 k\dy\d 2 k 2 dy 2 d 2 k\dy\d 2 k 2 dy 2 d 2 k\dy\d 2 k 2 dy 2 

This production cross section is the main formal result of this work. We analyze the properties of the cross section 
([39]) below. 

IV. LONG-RANGE RAPIDITY CORRELATIONS: AWAY-SIDE AND NEAR-SIDE; HBT 

CORRELATIONS 



Our goal now is to evaluate the correlations resulting from the two-gluon production cross section (|39|) . that is 
from the cross sections in Eqs. (|3"3")l and (f3"5)l . The first step is to evaluate the interaction with the target. We will 
be working in the quasi-classical MV/GM limit, where the interaction and, hence, the production cross sections (|33|) 
and (|38[) . are rapidity-independent. 

We begin with the cross section in Eq. (|33|) . Even in the quasi-classical limit, evaluation of the interaction in 
Eq. (|3"Tj) is calculationally intensive (though conceptually rather straightforward) . To simplify the calculation we will 
also employ the large- N c expansion: we assume that the nucleons in the nuclei are made out of an order- N 2 valence 
quarks (or gluons), such that the saturation scale, which in such case is proportional to Q 2 ~ ct 2 N 2 , is constant 
in the 't Hooft's large- N c limit. In the saturation physics framework such approximation was used in [63|] for the 
quadrupole operator giving a reasonably good approximation to the exact answer [66| . Similar approximations are 
frequently used (albeit, often implicitly) in applications of anti-de Sitter space/conformal field theory (AdS/CFT) 
correspondence to collisions of heavy ions modeled by shock waves [68l-l7l| . 

At the leading order in 1/N 2 expansion of Eq. (|3"3"|) the interaction with the target factorizes, as discussed around 
Eq. (|34|) . In addition, the cross section in Eq. (|38|) is 1/iV^-suppressed (as compared to the leading term in Eq. (|33p) 
and can be neglected at the leading order in N c . The correlation function is then given by Eq. ([T7| with the single 
gluon production cross section from Eq. $Z^. In the MV/GM approximation the gluon color dipole interaction with 
the target is [53| 



S G (xi,x 2 ,y = 0) = exp 



1 |_ |2 n 2 f x l+ x 2\ ,_ f 1 



l 1 * 1 -* 21 Q A-^P»K\x x -x 2 \K 



(40) 



with <3 s o the rapidity-independent gluon saturation scale in the quasi-classical limit evaluated at the dipole center-of- 
mass (x\ + x 2 )/2 and A an infrared (IR) cutoff. We see that the quasi-classical single gluon production cross section 



is rapidity-independent and, for unpolarizcd nuclei and for perturbatively large kr 54 1. is also independent of the 
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azimuthal angle <p °f the transverse momentum k of the outgoing gluon. We conclude that at the leading order in 
1 /N 2 in the quasi-classical approximation the geometric correlations are almost absent and the correlation function 
(flT]) is approximately zero. 

Non-trivial correlations can be obtained from Eq. Q33[) by expanding the interaction term to the first non-trivial 
order in 1 /N 2 . To this end we write a correlator of two Wilson line traces as 



1 



(n 2 iy 



(Tr[U Xl Ul 2 ]Tr[U X3 UlA) 



1 



(N 2 



(Tr[U Xl U X2 })(Tr[U X3 Uij) + A( Xl ,x 2l x 3 ,x 4 ) 



(41) 



where the correction to the factorized expression denoted by A, is order-l/A 2 . (Note that each adjoint trace in 
Eq. (|4Tj) is order-(A^ — 1), such that its left-hand side, as well as the first term on its right-hand side, are order-one 
in N c counting.) 

The leading order- 1/N 2 contribution to A is derived in the Appendix, with the result being 



A(x 1 ,x 2l x 3 ,x 4 ) 



(D 3 



N 2 



2e l 



2e L 



Di - D 2 (D 1 - D 2 ) 2 Dx - D 3 



(Dx - D 3 ) 2 



where we have defined 



Dx 
D 2 
D 3 



Q 2 



?s0 



35 1 - X 2 



35 1 - 35 3 



2 In 
2 In 



1 



Q 2 



\xx — X/x\ 2 In 



|35i - 35 2 |A 
1 

|35i - 35 3 |A 
1 

|35i — X/i\K 



(Dx - D 3 f 



\x 3 — x±\ 2 In 



\x 2 — 354I In 



1x2 — 35 3 | 2 In 



O 



1 



\x 3 — 354|A 



|35 2 - 35 4 |A 



|35 2 - 35 3 |A 



(42) 

(43a) 
(43b) 
(43c) 



assuming, for simplicity, that all the saturation scales are evaluated at the same impact parameter. 
Using Eq. (|41j) in Eq. (|33|) we see that the correlated part of the two-gluon production cross section is 



da. 



(corr) 
square 



a 2 C 



d 2 kxdyxd 2 k 2 dy 2 
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J d 2 Bd 2 bx d 2 b 2 Tx(B ~bx)Tx(B ~b 2 )d 2 xxd 2 yx d 2 x 2 d 



bx 



Vx - bx x 2 
\xx-bx\ 2 ' \yx-bx\ 2 \x 2 - 



V2 



1 2/2 



x [A(a5i,yi,35 2 ,y 2 ) - A(xx,yx,x 2 ,b 2 ) - A(xx,yx,b 2 ,y 2 ) - A(x l7 bx,x 2 ,y 2 ) - A(bx,yx,x 2 ,y 2 ) 
+A(xx,bx,x 2 , 6 2 ) + A(35i, 61, b 2l y 2 ) + A(b 1 ,y 1 ,x 2 , b 2 ) + A(bx,yx,b 2 , y 2 )] • (44) 

Eq. (|4^|) along with the expression for A in Eq. (|4"2")l arc our most complete results for the contribution to the two-gluon 
production cross section coming from Eq. (|33l) in the quasi-classical regime of the heavy-light ion collisions in the 
large- N c limit. The evaluation of the full Eq. (|44[) appears to be rather involved and is left for the future work. 

Instead we will expand Eq. (|42[) and use the result in Eq. (|44j) to obtain the correlated gluon production at the 
lowest non-trivial order. This result can be used to elucidate the structure of the long-range rapidity correlations, 
along with the comparison to the existing expressions in the literature. 

At the lowest non-trivial order in ZVs (that is, since each Di represents a two-gluon exchange with the target, at 
the lowest order in the number of gluon exchanges corresponding to kxr, k 2 T 3> Qso), Eq. (|42|) becomes 



A(35!, 35 2 , 353 , 35 4 ) 



(D 3 



2N? 



(45) 



Substituting this into Eq. (|44[) we can further simplify the expression by assuming that for the connected diagrams 
that contribute to the A's one has bx and b 2 pcrturbativcly close to each other, with the typical separation between 
these two impact parameters much smaller than the nucleon radius. Since the nuclear profile functions for a large 
nucleus do not vary much over perturbatively short distances we can put bx ~ b 2 rj b in the arguments of Tx's and 
Qso, where b = (bx + b 2 )/2. Defining A6 = bx — b 2 we can write d 2 bx d 2 b 2 = d 2 bd 2 Ab. The integral over A6 can be 
then carried out with the help of a Fourier transform 



1 



35I 2 In 



1 



MA 



r] 2 l 1 
2ir [ ' (P) 2 



(46) 
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used to replace (D3 — D 2 ) 2 in Eq. (|45|) and, hence, all A's in Eq. (|44|) as a double Fourier-integral: for instance 



A(x 1 ,y 1 ,x 2 ,y 2 ) 



d 2 ld 2 V 1 1 



2N 2 J (2tt) 2 (I 2 ) 2 (I' 2 ) 

x 

Qfo fd 2 ld 2 V 1 



e ii-(xi-x 2 ) 1 e i/-(yi-y2) _ e iJ-(asi-j/ 2 ) _ g i Z-(yi -sc 2 ) 

e -il' -(xi-x2) 1 g-ji'-(s/i— yz) _ e — tl'-(asi— #2) _ e — ii'-(yi— a) 3 ) 
1 



2iV 2 7 (2tt) 2 (Z 2 ) 2 (Z' 2 ) 2 



il-(x 1 -x 2 ) 1 gii-(yi-y2) _ gi/-(sci— 1/2) _ g»i'(j/i— x 2 ) 



-il' ■(x 1 -x 2 ) _|_ g— 2 l'-(yi -3/2) _ g— ii'-(sci— V2) _ g— ii'-(vi— S2) 

where we employed the substitution 

£ 1 = x 1 -6i, yi=yi-bi, x 2 = - 6 2) !/2=lte-&i- 



,i (l-l')-Ab 



(47) 



(48) 



Performing similar substitutions for all A's in Eq. (|4"4"]) . we can integrate over Afo, a?i, 0:2, j/i, 2/2, and Z'. After some 
algebra one arrives at (cf. (ill Il2|) 



(corr) 
square 



d 2 k\dy\d 2 k 2 dy 2 



LO 



d 2 



+ 



1 



(fel -l) 2 (k 2 + l) 2 (k 1 -l) 2 (k 2 -l) 



(49) 



where the subscript LO denotes the lowest-order cross section. 

Equation (|4TJ)) is illustrated in Fig. [5] by regular Feynman diagrams that contribute to its right-hand-side (cf. 
[Hill]): these dia grams are referred to as the "glasma" graphs in the literature. The momenta of the gluon lines are 
labeled in Fig. and the triple gluon vertices, marked by the dark circles, are the effective Lipatov vertices. These 
particular graphs contribute to the first and the second terms in the square brackets of Eq. (|49p correspondingly, and 
can be calculated by taking two Lipatov vertices squared. 



k 2 + 




k 2 -l 



O 
O 

o 
o 
o 
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FIG. 8. Examples of diagrams generating contributions to Eq. (|49[1 : the left panel represent the away-side correlations (the 
first term in Eq. Q49| l1. while the right panel contributes near-side correlations (the second term in Eq. (|49|l ). The f-channel 
gluon momenta flow toward the triple-gluon vertices to the left of the cut, and away from those vertices to the right of the cut. 



The obtained expression (|49| contains both the near-side and away-side azimuthal correlations [Hi, |12| : clearly the 
first term in the square brackets of Eq. (|49|) contains poles at I = k\ and I = —k 2 , which, after integration over Z, 
lead to a contribution^ 

~ WTW (50> 



2 Note that both terms in Eq. II49H also contain a pole at I = 0, which leads to a correlated contribution independent of the azimuthal 
angle between fci and kz- 
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characteristic of the away-side correlations. 

The second term in the square brackets of Eq. (|4"5)) has poles at I = fci and I = fc 2 , yielding a contribution 

1 



(fci - k 2 y 



(51) 



indicating near-side correlations. Note that all correlations are long-range in rapidity since the cross section 
rapidity- independent . 

Now we turn our attention to Eq. ([55)) . There the cross section contribution itself is l/iV 2 -supprcsscd as compared 
to the leading (uncorrelated) part of Eq. (|33f : hence we need to evaluate the interaction with the target in (|38j) using 
the large- N c limit. We work in the same quasi-classical MV/GM approximation, al ong with the large- N c limit. The 
fundamental (quark) quadrupolc amplitude was evaluated in this approximation in [63J (see Eq. (14) there), yielding 



'quark 



{x 1 ,x 2 ,x 3 ,x 4 ) = e Dl/2 



Do - Do 



e D 1 /2 _ p D 3 /2 



Dx-Dz 

Since the adjoint (gluon) quadrupolc (|36[) in the large- N c limit is simply 

Q{X 1 ,X 2 ,X 3 ,X 4 ) = [Q q uark(x 1 ,X2,X 3 ,X 4 )} 2 

we get 

D 3 - D 2 



Q(x 1 ,x 2 ,x 3 ,Xi) 



D, - D, 



p D 1 /2 _ e D 3 /2 



(52) 



(53) 



(54) 



Equations (|54j) and ()40|) , when used in Eq. (j38j , give us the remaining contribution to the two-gluon production cross 
section in the quasi-classical approximation, which has to be added to Eq. (|4"4"]) . 

Again the full expression (|3"5]l is hard to evaluate: instead, similar to Eq. (|4l?|) . we will consider the limit of large 
k\T = |fei| and k 2 T = \k 2 \: k\,k 2 ^> Q s o- Expanding Eq. (|5~4"]) in the powers of £Vs yields 

Q{x u x 2 , x 3 , x 4 ) = 1 + D x - D 2 + D 3 + 1 [2 D\ + D\ + 2 D\ - 3 D x D 2 - 3 D 2 D 3 + 3 D x D 3 ] + O (£>|) . (55) 

Using this result along with a similar expansion for Eq. f|40[) in Eq. p7p . and employing again the Fourier transform 
(|46p along with the substitution (j48j one can write 



Int crossed (xi,yi,bi,x 2: y 2 ,b 2 ) = Q% 
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(56) 



where the /<->!',!(-> — i, and Z' <-> —I' symmetries of the integrand were utilized to cast the expression in its present 
form. Using Eq. (|56p to replace the interaction with the target in Eq. Q38p and integrating over x x , x 2 , y\, y 2 , and, 
in some terms, over V , leads to the following result 
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where 
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(58) 
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and 
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HBT 



d 2 kidy\d 2 k 2 dy 2 



LO 
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(59) 



We will explain the origin of the "HBT" label on the cross section in Eq. ([59]) in a little while below: we defer the 
analysis of Eq. ([5U[) until then. 

First we concentrate on the expression (|58[) . We note that by shifting the integration momentum I one could reduce 
the second term in the curly brackets (together with the k 2 — ► —k 2 term) to that in Eq. (j4l?)) , thus doubling Eq. (|49l) . 
That term is still described by the diagrams of the type shown in Fig. [8] and gives the near- and away-side correlations 
in Eqs. {50) and ([ST]) . 
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FIG. 9. An example of a diagram giving a contribution to the first term in the curly brackets of Eq. ([58 



The first term in the curly brackets of Eq. (|58[) corresponds to a different class of diagrams, one of which is shown 
in Fig. [9] One can see that the diagram in Fig. |9] contributes two non-forward "squares" of the effective Lipatov 
vertices. An analysis of the poles of the integral in the first term of Eq. (|58p shows that it contains similar near- and 
away-side correlations to that in Eqs. (|50[) and (|51[) . though with an additional enhancement due to a prcfactor, 



[2(fc 1 -fc 2 ) 2 -fc 2 fe2 
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(fci-fca) 2 (fei+fe 2 ) 2 



(60) 



The origin of this contribution to correlations is in the non-forward squares of the Lipatov vertices, which is similar to 
the mechanism for generating long-range rapidity correlations proposed in [39j using two BFKL ladders with non-zero 
momentum transfer. Note, however, that we do not obtain the correlation proportional to the first power of k\ ■ k 2 
advocated in [39(, perhaps due to the lowest-order nature of the result ([55)) . It appears that the correlations ([6TJ]) 
resulting from the lowest-order diagrams like the one depicted in Fig.[9]have not been considered in the literature yet. 

Note that the away- and the near-side correlations enter Eqs. ([4"9")) and ([55j) on equal footing: in fact, one could 
be obtained from another by a simple k 2 —> —k 2 substitution in cither of those expressions. While evaluating the 
integrals over I in Eqs. (|49p and ([58[) analytically appears to be rather algebra-intensive, we do not need to do this 
to observe that, once the /-integral is carried out for one of the terms in the square brackets of Eq. (j49[) , the answer 
for the other term is obtained by substituting k 2 — > —k 2 in the result. Equation ([58[1 simply contains an additive 
k 2 — > — k 2 term. 

Since the correlated cross sections ([491) and (|58[l are sums of two terms related by the k 2 — > — k 2 substitution 
and are symmetric under the k\ f-> k 2 interchange, we conclude that they are functions only of even powers of k\ 
and k 2 , that is functions of kf, k 2 , (ki ■ k 2 ) 2 , and possibly (fei x k 2 ) 2 . Clearly this implies that the Fourier series 
representation of Eqs. (j4U)) and ([55)) would only contain even cosine harmonics of the azimuthal angle, that is 
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where A0 — (f>i ~ 4> 2 is the angle between momenta fci and fc 2 while kyr = |fei| and fc 2 T = I&2I- Here c„(fci,A; 2 ) are 
some coefficient to be determined by an exact calculation. It is quite interesting that only even harmonics contribute 
to the correlated cross section in Eq. (J6TJ . Let us stress here that we have not made any assumptions about the 
centrality of the collision: we do not have an almond-shaped overlap of the two nuclei. In fact we integrate over all 
impact parameters B. (Also the impact parameter dependence factorizes from the rest of the expression.) Hence the 
correlation in Eq. (|61[) is not caused by the geometry of the collision. 

To construct the correlation function one has to use Eqs. (|4"!)1) and (f5"8"l) in Eq. (JT5J) . In the latter, the uncorrelated 
two-gluon production would dominate in the denominator of the normalization factor, that is, in the denominator of 
the first factor on its right-hand side. The singlc-gluon production cross-section ([29]) at the lowest order is equal to 
(see [l6[ and references therein) 
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Using this along with the sum of Eqs. (14TH) and (|5"8"|) in Eq. (|T51) and taking the large- N c limit in Eq. (|62[) we get 
1 Jd 2 Bd 2 b[T 1 (B-b)} 2 Q%(b) 



(62) 
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(63) 



To evaluate this further let us assume for a moment that both colliding nuclei are simply cylinders oriented along 
the collision axis, such that the nuclear profile functions of the nuclei are Tj(6) = 2 pRi 8(Ri — b) where i — 1,2 labels 
the nuclei, p is the (constant) nucleon number density in the nucleus, i?i and R 2 are the radii of the projectile and 
target nuclei, and the cylindrical nucleus is assumed to have length 2Ri along its axis. Assuming that both nuclei are 
large and neglecting the edge effects, we can neglect the bi and b 2 dependence in these nuclear profile functions of 
both nuclei. Since the gluon saturation scale in the MV model is Q 2 = 4ira 2 T 2 (b) we obtain in the R\ <C R 2 limit 
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(64) 



Indeed the correlator is suppressed by a power of N 2 and a power of the cross sectional area of the projectile nucleus 
ttR 2 , as discussed in the literature plTfl4l l40j . 

While our conclusion about the correlator (|64|) and the cross section Eq. ([61) contributing only to even Fourier 
harmonics has been verified above at the lowest order only, it is true for the full two-gluon production cross section in 
heavy- light ion collisions (|39|) . This can be seen by noticing that fc 2 — s- — fe 2 substitution in Eq. (|44|) does not change 
the cross section, as it is equivalent to the x 2 o y 2 interchange of the integration variables. The integrand of Eq. (|44|) 
(or, equivalently, of Eq. (|33p ) is invariant under x 2 <-> y 2 interchange since the gluon is its own anti-particle such that 
Tr [U x C/JJ = Tr \U y U£\ . The expression in Eq. (|58|) is explicitly invariant under fe 2 — > —k 2 substitution. Hence the 
net correlated cross section (|39|) is an even function of fc 2 . Note also that Eqs. (|33|) and (|38|) are fei fe 2 -symmetric: in 
Eq. (|33p and in the term arising from the first exponential in Eq. (|38j) the symmetry is a consequence of the symmetry 
of the integrand under the simultaneous X\ •H- x 2 , y\ O y 2 , and b\ o b 2 interchanges. The term multiplying the 
second exponential in Eq. (|38[) is fei -f-> fc 2 -symmetric due to the X\ y\, x 2 -H- y 2 symmetry of the integrand, which 
follows from the following property of the quadrupole operator: Tr [U Xl U X2 U^ ] = Tr [U V2 Ul 2 U yi U Xl ]. We 
see that the net cross section (|39p is decomposable into a Fourier series with even harmonics only. Therefore, the 
correlation function in the heavy-light ion collision can be also written as an even-harmonics series 
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with some coefficients d n . This conciusion also seems to hold in the case of classical gluon fields produced in a collision 
of two heavy ions, as can be seen from the result of the full numerical simulation of two-gluon production in heavy 
ion collisions due to classical gluon field carried out in [72j]. The correlators in Figs. 9 and 10 of [72j do appear to have 
similar- looking maxima at A</> = and A0 = 7r (within the accuracy of the numerical error bars), though a more 
careful analysis is needed to figure out if our conclusion is true for two colliding heavy ions. 

To better visualize the correlator let us envision a toy model in which multiple rescatterings (and other saturation 
effects) regulate the singularities at k\ = ±k 2 by the saturation scale Q s q in such a way that the correlator can be 
modeled as proportional to 

Ctoy m0 del(fci^^2) ~ {ki k ] )2 - + {ki+k ] ? + Qlo (66) 

This correlator is plotted in Fig. [TU] using arbitrary units along the vertical axis as a function of the azimuthal angle 
A0 between k\ and k 2 for k\ = k 2 = Q s o- The shape illustrates what the full correlation function may look like, 
having identical near- and away-side correlation peaks. 



am 




FIG. 10. A toy azimuthal two-gluon correlation function motivated by the calculation in this Section. Vertical scale is arbitrary. 



Our conclusion in this Section is that the saturation/CGC dynamics in nuclear collisions appears to generate the 
long-range rapidity correlations which have identical maxima at both A<j) = and A(f> = n in the azimuthal angle. 
Such correlations are non-flow in nature, since they do not arise due to almond-shaped geometry of the collisions. (In 
fact the correlations should persist for the most central heavy ion collisions, though they should be suppressed by an 
inverse power of the overlap area as discussed above and in [llT - fTil |40||.) However, since the elliptic flow observable v 2 
(even in the reaction-plane method) is determined from two-particle correlations, it is possible that the correlations 
discussed here may contribute to the elliptic flow (and to higher-order even- harmonics flow observables V2n) measured 
experimentally. It is, therefore, very important to experimentally separate our initial-state correlations from the 
late-time QGP effects. 

Naturally the cumulant analysis [73|, [74[ is likely to remove these two-particle non-flow correlations from the flow 
observables. However, the effectiveness of the cumulant method may again depend on the collision geometry. Let us 
illustrate this by considering the fourth order cumulant for elliptic flow, defined by [73|, [zH 

c 2 {4} = (e 21 (^i+*2-^-04)^ _ ^ e 2i(0i-<fe)\ ^ e 2i(0 2 -0 4 )^ _ ^ e 2i(0!-0 4 )^ / e *i{<h-<t>a)\ (67) 

where the angle brackets denote event averages along with the averaging over all angles <f>i , <p2 , 03 , and </>4 of the four 
particles employed in the definition. Using the lowest-order correlations employed in arriving at Eq. (|63[) one can 
straightforwardly show that the cumulant due to these correlations only is proportional to 

m{ Jd 2 Bd 2 b 1 d 2 b 2 [T\{B bx)] 2 [Ti(B- b 2 )} 2 Q^jbx) Q%(b 2 ) 

C2t j liO J d 2 Bd 2 bl d 2 b2 d 2 h d 2 h Tl{B _ bl) Tl[B _ b2 ) Tl{B _ b3) Tl( _B _ b4 ) Q2 o(bi) g2 o(b2) Q^ Q (b 3 ) Q 2 (b 4 ) 

J d 2 B d 2 b [Ti(B - ft)] 2 Qs (b) ] 2 



J d^B d 2 &! d?b 2 T X {B ~b 1 )T 1 (B^b 2 ) Q 2 (bi) Q 2 (f> 2 ) 

if the 4-particle correlator in Eq. (|67|) (the correlated part of the first term on its right-hand-side) is due to the 
pairwise correlations of particles 1,3 and 2,4 or 1,4 and 2,3. Clearly, in the general case, C2{4}| io from Eq. ([55)1 
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is non-zcro0 Moreover, the fourth order cumulant would contain the azimuthal angles dependence resulting from 
Eq. (j6"3")l but not shown explicitly in Eq. We see that the geometric correlations from Sec. HT1 may prevent 

complete removal of these non-flow correlations in the cumulants. Just like with the correlator of Eq. ([20]) . the above 
non- flow correlations can be completely removed from the cumulants for the fixed impact parameter B: if we fix B in 
Eq. (|68|) (that is, remove all the d 2 B integrations from it), we get C2{4}| LO = 0, which is exactly what the cumulant 
is designed to do [lH, [74| — completely cancel for non-flow correlations. Note that since, in the actual experimental 
analyses one effectively integrates over B in a given centrality bin, it is possible that some non-flow correlations (|63j) 
would remain in the cumulant (|68j) . Even fixing \B\ precisely and integrating over the angles of B may generate a 
non-zero contribution of these correlations to the cumulant. The question of the interplay of the true QGP flow and 
the non-flow correlations discussed here has to be resolved by a more detailed numerical study. 

To clarify the physical meaning of the cross section obtained in Eq. (|59|) we again consider a collisions of cylin- 
drical nuclei. With this simplification we can consider one of the terms in the interaction with the target, say the 
Q(xi, yi, x 2 , yz) in Eq. (|37p . Among many terms which contribute to the quadrupole amplitude in Eq. (|55[) . there is 
a term 



Q(xi,yi, x 2 ,y 2 ) 
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(69) 



Using Eq. (|69|) in Eq. (|38|) and changing the interaction variables to those defined in Eq. (|48|) we see that the remaining 
b\ and b 2 integrals, after shifting those variables by B, become simple Fourier transforms of the projectile nucleus 
profile function Ti(6), 
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(70) 



In the limit of a sufficiently large projectile nucleus the Fourier transforms give a delta-function, yielding 

d<7HBT 



d?kidyid?k 2 dy 2 



5\kx-k 2 ) + 5\kx + k2) 



(71) 



in agreement with the result in Eq. (|59[) . Our present estimate shows that a more careful evaluation of Eq. (|70p would 
give a smother peak in \k\ — k 2 \ at k\ = k 2 with the width determined by the inverse radius of the projectile nucleus 
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FIG. 11. An example of a lowest-order diagram generating HBT-type correlations. 



The first term on the right-hand-side of Eq. ((71]) (or Eq. (|70l) ) has the trademark form of the Hanbury Brown-Twiss 
(HBT) correlations [EH , which are widely studied in heavy ion physics [75l - l77| . These correlations are normally local 
both in rapidity (y\ = y 2 ) and in the transverse momentum k\ = k 2 : however, in Eq. (|59j) (or, in Eq. (|7ip ) we only 



3 Note, however, that for the heavy-light nuclear collision case considered here, A\ <JC Ai, and for cylindrical nuclei, one gets c 2{4}| i0 = 0. 
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have locality in transverse momenta. This can be explained by the fact that the extent of the interaction region in 
the longitudinal direction, commonly labeled Ri ong , is very small in our case due to the extreme Lorentz-contraction 
of the two colliding nuclei leading to smearing of the longitudinal HBT correlations. 

We conclude that the correlations resulting from Eq. (|38|) include HBT. Diagrammatic representation of the HBT 
correlations was discussed earlier in [78j . An example of the diagram giving rise to the correlations in Eq. (|59j) is 
shown in Fig. ITT1 

Interestingly though, due to the f-> — hi symmetry of the two-gluon production cross section (|38p . the HBT peak 
at fci = &2 in Eq. ([59]) is accompanied by an identical peak at k\ = — k^ resulting from S 2 (ki + A^). We thus obtain 
a back-to-back HBT correlation resulting from multiple rescatterings in nuclei. (The origin of our back-to-back HBT 
correlations is different from that of the back-to-back HBT-like correlations proposed in [79l |80||.) 

Note that it is possible that the process of hadronization would affect the phases of the produced gluons, possibly 
destroying the perturbative HBT correlations of Eq. (|38p and replacing them with the HBT correlations of the 
non-perturbative origin. The same hadronization process may also destroy the back-to-back HBT correlations from 
Eq. dMD- 



V. CONCLUSION 

Above we have derived the cross section for two-gluon production in the heavy-light ion collisions. The cross section 
is given by Eq. (|39[) . with the two contributions shown in Eqs. (|33[) and (|38[) . Concentrating on gluon production with 
a long-range separation in rapidity of the two gluons we evaluated the interaction with the target in Eq. (|33[) using 
the MV/GM approximation in Eq. (jl2"j) . The interaction with the target in Eq. in the MV-model framework can 
be obtained from Eqs. (gDJ) and 




FIG. 12. Two-gluon production in a standard pQCD formalism, with the diagram on the left generating back-to-back correla- 
tions, and the diagram on the right responsible for the near-side collinear correlations. 

Analyzing the corresponding correlation function we identified four main types of correlations: (i) geometric corre- 
lations, (ii) HBT correlations, (hi) away-side correlations, and (iv) the near-side correlations. HBT correlations (ii) are 
local in rapidity and thus are outside the main scope of our work here: however, we find that HBT correlations arising 
in our calculation are accompanied by a back-to-back HBT peak, as follows from Eq. (f5T)|) . Geometric correlations 
(i) arise in Eq. (|17[) due to nucleons being confined within the same nucleus: more work is needed to see how strong 
they are and what role they are playing in experimental data. 

Concentrating on the long-range rapidity correlations (iii) and (iv) we have shown that the two correlations have 
identical azimuthal profiles centered around A<fi = tt and A4> = correspondingly. This result is true to all orders 
of the calculation in the heavy-light ion collision considered here. Let us stress here the difference between this 
result and the regular jet correlations: indeed a two-hadron production cross section calculated in perturbative QCD 
(pQCD), illustrated by the examples in Fig. [TSJ contains both the back-to- back and near-side correlations. However, 
while the jet back-to-back correlations are long-range in rapidity, the jet near-side correlations are local in rapidity, 
being due to collinear jet showers around the trigger hadron (see (8ll . |82| for an analysis of mini-jet correlations in 
heavy ion collisions). Hence in standard pQCD formalism the strengths along with rapidity and azimuthal shapes of 
the back-to-back and near-side correlations are different. This is in stark contrast to our correlations (iii) and (iv), 
which arc both long-range in rapidity, and have identical amplitudes and rapidity /azimuthal shapes. The correlations 
(iii) and (iv) arise due to multiple rescatterings in both nuclei and are, therefore, typical of the ion-ion collisions, 
weakening in the limit when one (or both) of the ions is (are) small (a proton). It is possible that this correlation 
would complicate the experimental extraction of the elliptic flow observable V2 due to the QGP flow in heavy ion 
collisions, though it may also be that the correlation is too weak to present a significant background. More detailed 
numerical analysis of our result is needed to clarify this point. 
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VI. APPENDIX 



Our goal here is to calculate the following object 

(jy2 \ 1)2 (Tr[U ai uU Tr[U X3 Utj) (Al) 

in the quasi-classical MV/GM approximation at the lowest non-trivial order in 1/N% expansion in the 't Hooft's 
large- A c limit. We assume that the saturation scale is A c -independent, that is Q 2 ~ (N c )°, due to an order- N 2 
number of "valence" partons in each nucleon in the target nucleus. 




FIG. 13. Diagrams contributing to the correlator in Eq. (|A1|I in the quasi-classical approximation. Vertical dashed lines denote 
t-channel gluons, while the shaded ovals represent the nucleons in the target nucleus. 



A sample of the diagrams contributing to Eq. (|A1[) is shown in Fig. Q21 where the t-channel gluons are denoted 
by dashed lines to simplify the picture. The calculation is straightforward [Hj: one simply has to exponentiate the 
two-gluon exchange interaction with a single nucleon (the nucleons are denoted by shaded ovals in Fig. I13|) . The only 
complication is that, unlike the dipole amplitude calculated in (53|, the interaction now, for the double-trace operator 
(jAlj) . is a matrix in the color space. A double trace operator like (|Alj) with the Wilson lines in the fundamental 
representation was calculated earlier in [85| (for similar calculations see also [65l IoTI |86|). 

To exponentiate the matrix we have to choose a basis in the color space of four s-channcl gluons in Fig. 1131 clearly 
the net color of the four gluons is always zero. The color states of the four gluons can be classified according to the 
color states of the top two s-channel gluons, since the color of the bottom pair of s-channcl gluons is determined by 
requiring net color-neutrality of the four s-channel gluon system. The colors of a gluon pair can be decomposed in 
the following irreducible representations 

(n 2 - 1) g> (n 2 - 1) = Vi ® v 2 e v 3 e v A e v 5 e v 6 e v 7 

= i e (n 2 - D e N ^ - » + 1} e N2 ^ + ^ - 1} e (ni - 1) e - m - 4 ) ffi w - m - 4 ) . 

4 (A2) 

Here we are following the notations introduced in [87|, see page 120 there. Labeling the colors of the four s-channel 
gluons in an arbitrary color state by a, b, c, and d as shown in Fig. [T3] we define the color states corresponding to 
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representations V±, Vi, and V5 by 
1 



Pi = 



m - 1 



|P 2 ) = 



d abe d cde , |P 5 ) = 



be rede 



y/N* - 1 



jabe £ 



(A3) 



where we differ from Pj's in [87j by prcfactors, since here we demanded that our color states are normalized to one, 



(Pi\Pj) = 5 V} . Other color states can be constructed as well [87J, but we will only need the states in Eq. (|A3| for the 
calculation below. 

We denote by M the interaction with a single nucleon by a two-gluon exchange: it is a matrix in the color space 
of the s-channel gluons. Since in the correlator (|Alj) the top (bottom) two s-channel gluons are in the color-singlet 
state both before and after the interaction, we write 



1 



(Tr[U Xl Ul 2 ]Tr[U X3 Ul}) 



(Pi\e M \P 



1 ■ 



(iV| - l) 2 

Expanding the exponential in a power series and inserting unit operators 1 = \Pi] 

1 



' 11 



(A4) 

(Pi I between all the M's yields 

(A5) 

(A6) 

All one has to do now is to find the matrix M from Eq. (|A6p . exponentiate it, picking up the "11" matrix element of 
the exponential. Such a calculation, while straightforward, is rather involved: here we will utilize the large- N c limit 
to construct an approximate result. 

Calculating some of the elements of the matrix M and evaluating the A^-order of the remaining matrix elements 
yields 



(iV c 2 - I) 2 

where the 7x7 matrix M is defined by its elements, 

Mij = (Pi\M\Pj) 



( A 



M = 







D3-D2 








\D 1 + \(D 2 +D a 
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with Di defined in Eqs. (|43j) and ellipsis denoting the matrix elements we do not need to calculate as they are at most 
order-1 in N c counting. From Eq. (|A7[) we sec that if we start in the color-singlet state \P±) for the top (bottom) gluon 
pair, a single interaction can either leave the two gluons in the color-singlet state, or can flip them into a color-octet 
state I-P5). The latter transition comes in with an order-l/A^ c suppression factor. In order to evaluate (|A1|) we have 
to start and finish with a color-singlet state: to order-l/A", 2 we may have at most two such transitions: \Pi) — > IP5) 
and the inverse, |Ps) — > |Pi). Once the system is in the color-octet IP5) state, it can continue its random walk 
through color space: however, if we want to keep the calculation at the ordcr-l/A^ 2 , the interactions between the two 
transitions should be leading-order in N c . This is why only the leading- A^-order matrix elements M22, M25 = -M52, 
and M55 contribute in the 6x6 matrix with i,j = 2,...,7: we do not need to calculate the l/A^-suppressed 
elements in Eq. (|A7|) or the elements denoted by ellipsis which can not contribute. 

Exponentiating the matrix M from Eq. (|A7j) . picking up the "11" element of the obtained matrix and expanding 
the result to order- yields 
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(iV 2 - 1)^ 
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Finally, noting that (see Eq. (|gQjl) 
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(N*-iy 



(Tr[U Xl Ul \) (Tr[U X3 Ul }) = e Dl 



(A8) 



(A9) 
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and using Eq. (|41~j) we obtain Eq. (|42|) . 
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